Relative proton and 7 widths of astrophysically important states in 30 S studied in the 

/?-decay of 31 Ar. 
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Resonances just above the proton threshold in 30 S affect the 29 P(p, 7) 30 S reaction under astro- 
physical conditions. The (p, 7)-reaction rate is currently determined indirectly and depends on the 
properties of the relevant resonances. We present here a method for finding the ratio between the 
proton and 7 partial widths of resonances in 30 S. The widths are determined from the /32p- and 
/3p7-decay of 31 Ar, which is produced at the ISOLDE facility at the European research organization 
CERN. Experimental limits on the ratio between the proton and 7 partial widths for astrophysical 
relevant levels in 30 S have been found for the first time. A level at 4689.2(24) keV is identified in the 
7 spectrum, and an upper limit on the Fp/F-, ratio of 0.26 (95% C.L.) is found. In the two-proton 
spectrum two levels at 5227(3) keV and 5847(4) keV are identified. These levels were previously seen 
to 7 decay and upper limits on the F-,/F p ratio of 0.5 and 9, respectively, (95% C.L.) are found, 
where the latter differs from previous calculations. 

PACS numbers: 26.30.-k, 23.40.Hc, 27.30.+t 
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I. INTRODUCTION 

Detailed knowledge of the energy levels of exotic nuclei, 
especially the ones just above the thresholds for particle 
emission, is important for understanding astrophysical 
processes such as explosive hydrogen burning. 30 S is sit- 
uated close to the proton drip line and is produced in the 
29 P(p, 7) 30 S reaction in the rp- and op-process in type I 
x-ray bursts [H HJ. The relatively long life time of 30 S 
makes it a critical waiting-point nucleus for these pro- 
cesses 0. 

The 29 P(p, 7) 30 S reaction is also interesting for the 
study of presolar dust grains. The most extensively stud- 
ied grains are SiC grains, because they are relatively 
abundant. A small fraction of these have been suggested 
to originate from classical novae [H . They are character- 
ized by low 12 C/ 13 C and 14 N/ 15 N ratios, high 30 Si/ 28 Si 
ratios, and 29 Si/ 28 Si ratios close to or lower than terres- 
trial values. The silicon isotopic abundance can provide 
information on the dominant nuclear synthesis paths fol- 
lowed by the thermonuclear runaway, which sets in near 
the base of the accreted layers from a main sequence star 
onto a white dwarf in a binary system [5|. In order to 
understand the origin of the isotopic ratios observed, the 
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processes that create and destroy the different silicon iso- 
topes have to be well understood. One of these is the 
29 P(p, 7) 30 S reaction. If this reaction is faster than the 
f3 + decay of 29 P, the amount of 30 Si would increase and 
the amount of 29 Si would decrease 0. 

The Gamow window of the 29 P(p, 7) 30 S reaction for 
temperatures relevant for astrophysics spans 100 keV to 
HOOkeV. This, along with the proton separation energy 
of 30 S being 4395.6(7) keV 0], implies that the levels in 
30 S interesting for astrophysics lie below 6 MeV. 

Iliadis et al. Q predicted that the reaction rate of the 
29 P(p, 7) 30 S reaction was dominated by two resonances 
in 30 S with spins 3 + and 2 + and excitation energies of 
4733(40) keV and 4888(40) keV, respectively, which had 
not been observed at that time. The first experimen- 
tal evidence came through studies of 32 S(p, i) 30 S by Bar- 
dayan et al. [9| and Setoodehnia et al. flfl. The most 
recent studies of 30 S can be found in Rcfs. while 
the results from previous experiments are combined in 
Ref. pj|. 

The reaction rate is calculated using the energy, spin, 
and the proton and 7 partial width of the relevant 
resonances. The improvements on the energies of the 
3 + and 2 + levels have reduced the uncertainties of the 
29 P(p, 7) 30 S reaction rate significantly for the relevant 
temperatures. The calculations made by Setoodehnia et 
al. [I4I shows that the uncertainties are now so small that 
they no longer significantly influence the silicon abun- 
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FIG. 1. (Color online) The p decay of 31 Ar. Different proton 
decays are drawn as an illustration. The marked proton lev- 
els in 30 S correspond to Fig. [3] and the marked 7 transitions 
correspond to Fig. [4] 



dances. These calculations, however, use proton and 7 
partial widths, which are calculated based on the shell 
model and comparisons with the mirror nucleus. Exper- 
imental values for these would clearly be preferred. We 
present here a method for finding the ratio between the 
proton and 7 partial widths for the resonances of astro- 
physical interest using the decay of 31 Ar. Until now the 
preferred method for studying these resonances has been 
by the use of reaction experiments. The results presented 
here open up a new approach. Further results from our 
experiment will be published separately. 

The experiment is described in SecHJ Sec. lIIII presents 
the results of the analysis of the low-lying states of 30 S 
and compares then with results in recent papers. It in- 
cludes the method for finding the ratio between the pro- 
ton and 7 partial widths. Finally, Sec. IIVI summarizes 
the main results. 



II. THE EXPERIMENT 

In this experiment, the resonances in 30 S are studied 
via the /3-delayed proton-7 decay and the /3-delayed two- 
proton decay of 31 Ar. The /3-delayed two-proton decay 
is known to be mainly sequential [ill ]. A partial decay 
scheme of 31 Ar is shown in Fig. [TJ The experiment was 




FIG. 2. (Color online) The experimental setup used for the 
experiment. The beam enters between DSSSD 5 and 6 and 
is stopped in a foil mounted on a small metal holder entering 
between DSSSD 3 and 5. Two clustered germanium detectors 
situated outside the cube behind DSSSD 3 and 4 are shown 
on the upper drawing. The top of the cube with three of the 
DSSSDs is lifted, following the dotted black line, for better 
visualization on the bottom drawing. 



optimized for detection of the delayed two-proton decay 
through a compact setup including detectors capable of 
stopping high-energy protons. This gave substantial (3 
background at low energy, which along with electronic 
noise made it challenging to identify low-energy protons. 

The radioactive 60-keV 31 Ar beam used in the exper- 
iment was produced at the ISOLDE facility at the Eu- 
ropean research organization CERN using a CaO target 
and a versatile arc discharge plasma ion source 11711 cou- 
pled to the General Purpose Separator (GPS) [18]. The 
beam is produced by irradiating the target with short 
high- intensity proton pulses at low repetition rate [l8j]. 
An average 31 Ar yield of about 1 per second was obtained 
for a runtime of about 7 days. 

The 60-keV beam was collected in a 50-^tg/cm 2 car- 
bon foil situated in the middle of the detector setup 
consisting of the silicon cube detector [ll|, containing 
six double sided silicon strip detectors (DSSSDs): one 
69 /jm (1), one 494 /im (5), and four close to 300 fim (2- 
4, 6). The setup can be seen in Fig. H The DSSSDs 
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are segmented into 16 strips in the front and in the back 
side, each 3 mm wide and 0.1 mm apart. Behind the thin 
DSSSD (1) and three of the four 300-/«n DSSSDs (2, 
3, 6) 50 mm x 50 mm unsegmented silicon pad detectors 
(thickness close to 1500 /Ltm) were placed, which enables 
particle identification. In some part of the analysis only 
the four DSSSDs backed by a pad detector are used. For 
detection of 7 rays, two cluster detectors, each consisting 
of three germanium crystals, from MINIBALL [2(| were 
placed outside the chamber containing the silicon cube. 

The energy and geometry calibration of the DSSSDs 
are made with 33 Ar produced from the same target-ion 
source unit as 31 Ar. A number of runs were made with 
33 Ar during the experiment to monitor if the energies 
drifted. 33 Ar decays by /3-delayed proton emission with 
well-known proton energies. For calibration of DSSSD 3- 
6 three known levels of 33 C1 were used: 3971.24(19) keV, 
4112.34(20) keV, and 5548.6(5) keV 0, corresponding 
to proton energies of 1642.8(5) keV, 1779.5(5) keV, and 
3171.9(7) keV. DSSSD 2 has several broken strips and 
thus fewer counts. The protons from the two lowest 33 C1 
resonances used could not be identified here. Instead, 
two higher-lying proton peaks are used. Their energies 
were found in the calibrated proton spectra of DSSSD 
3 and 6, which had a better energy resolution, giving 
proton energies of 2479.2(20) keV and 3856.9(20) keV. 
Due to the thickness of DSSSD 1 protons with energies 
above 2.50 MeV were not stopped inside the detector. 
For calibration of this detector the 33 C1 resonances at 
3971.24(19) keV and 4112.34(20) keV were used together 
with two proton peaks found from DSSSD 3 and 6 at 
1320.1(20) keV and 2479.2(20) keV. 

The solid angle of the DSSSDs is 43 % of 4tt, or, if the 
two detectors without backing are disregarded, 27% of 
4tt. Protons with energies below 500 keV can be stopped 
inside the collection foil depending on their emitted an- 
gle. The low-energy protons, which are interesting for 
this work, have an energy around 280 keV. The angular 
coverage of these protons is 20% of 4ir, when the two 
detectors without backing are disregarded. 

The half-life of 31 Ar is just 14.1(7) ms [lj]. The beam 
gate for 31 Ar was therefore only open for 100 ms after 
a proton pulse has hit the target. In addition to this a 
software time window from 5 ms to 100 ms after proton 
impact was used. 

The pad detectors are calibrated with a 148 Gd source 
and a triple a source consisting of 241 Am, 239 Pu, and 
244 Cm. 

The Ge detectors were calibrated using first a 137 Cs 
and a 60 Co source and then improved using a 152 Eu 
source, and lines from the decays of 32 > 33 Ar and 16 > 18 N 
were recorded online. The absolute 7 efficiency was found 
using a relative efficiency curve determined in a slightly 
different detector configuration (using four different 7 
sources: 152 Eu, 60 Co, 207 Bi and n Be) and an absolute 
measurement with an 152 Eu source. The result, using 



M 

o 



10 



o 
U 




10 



Excitation energy in S (MeV) 

FIG. 3. (Color online) The excitation energy of 30 S calculated 
from the energies of the two protons from the sequential two- 
proton decay as described in Ref. [l(|. The energies of the 
marked peaks can be found in Table U 
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FIG. 4. (Color online) The summed 7 spectrum over all crys- 
tals gated on protons from the decay of 31 Ar. The numbers 
correspond to transitions in 30 S and the Roman numbers cor- 
respond to transitions in 29 P. The energies of the marked 
peaks and their relative intensities can be found in Table [III 



the formula in Ref. , is 
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with an estimated uncertainty of 10 %. 

III. RESULTS AND DISCUSSION 

Table fl] summarizes the latest published results on the 
low- lying levels of 30 S. Included are also the energies 
measured in this experiment. Fig. [3] shows the levels 
above 5 MeV found, in the present work, using the en- 
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TABLE I. Energy and spin of the S levels below 6MeV for the present work and recent previous work. The proton separation 
energy of 30 S is 4395.6(7) keV 0. 



Setoodehnia et al. [12] 




Lotay et al. 


[13] 


Almaraz- 


-Calderon et al. 1 


'1 A~\ 

14 


Present work 


28 Si( 3 He,n 7 ) 30 S 


32 S( P ,t) 30 S 




28 Si( 3 He,n 7 ) 30 S 


32 S(p,t) 30 S 


28 Si( 3 He,n) 30 S 




31 Ar(/3+)(p) 30 S 




P flrpVl 




E x (keV) 




E x (keV) 


E x (keV) 




(keV) 


g.s. 




0+ 


g.s. 


0+ 


g.C5. 


g.b. 




g.S>. 


2210.6(3) 


2208(3) 


2+ 


2210.1(1) 


2+ 


2208 5(221 


2200(21 01 




2210 2C111 


3403.6(6) 




2+ 


3404.1(1) 


2+ 


3405 S?1 21 






3404 4(1 61 


3667.0(5) 




0+ 


3668.0(4) 


0+ 


3677.3(70) 


3600 ('2601 






Q(^7(^ C\( A \ 
oD (K).\)\Ql ) 


oDol^o ) 


1+ 


3677.1(4) 


1+ 










4688.1(4) 


4688(2) 


3+ 


4687.6(2) 


3 + 


4682.5(57) 






4689.2(24) 


4809.8(5) 


4812(2) 


2+ 


4808.7(3) 


2+ 










5132.3(5) 




(4+) 


5132.1(1) 


4+ 


5130.0(18) 




(4+) 






5225(2) 


(0+) 


5218.8(3) 


3+ 


5217.8(28) 


5200(44) 


(0+) 


5227(3) 




5315(2) 


(3-) 




(3-) 


5312.1(20) 




(3-) 






5393(2) 


3+ 




(2+) 


5382.0(7) 


5400(43) 


(2+) 


5392(4) 




5849(2) 


(2+) 


5848.0(4) 


4+ 


5835.5(13) 




(4+) 


5847(4) 




5947(3) 


(4+) 















TABLE II. The energies and relative intensities of the 7 lines 
in Fig. [4] normalized to the transition from the first excited 
state to the ground state of 30 S. 



Peak 


Energy (keV) 


Intensity 


1 


1194.2(11) 


35(4) 


2 


2210.1(11) 


100(8) 


3 


2478.9(21) 


16(4) 


4 


3407(7) 


22(7) 


I 


1383.7(14) 


19(4) 


II 


1957(4) 


10(4) 



ergies of the two protons from the sequential two-proton 
decay as described in Ref. [la ]. The two levels above 
6 MeV, seen in Fig. 02 corresponds to previously identi- 
fied levels [l6[ and will be treated in an other publica- 
tion. Due to contamination from /3 particles and elec- 
tronic noise, the levels below 5 MeV cannot be identified 
in the two-proton decay. These levels are known to de- 
cay mainly by 7 emission and are thus identified from the 
proton-gated 7 spectrum presented in Fig. S) The rela- 
tive intensities of the 7 lines are given in Table [TTJ Only 
the 1194.2- and the 2210.1-keV lines have previously been 
observed in the decay of 31 Ar [23[. The 1194.2- and the 
3407-keV lines both correspond to decays of the second 
excited state in 30 S. Their relative intensity is 0.63(21), 
which is consistent within 2ct with the value of 0.25(4) 
measured in Ref. [24j. However, the difference of these 
values might imply a double structure of peak 4 in our 
spectrum. In the following we focus on the relative 7 and 
proton decay widths of the 4689.2-, 5227- and 5847-keV 
levels. An overview of previously calculated r 7 /r p ratios 
can be found in Table UTTl 



The 4689-keV level 

The level at 4689.2(24) keV is observed for the first 
time in the decay of 31 Ar. It is identified in the 7 spec- 
trum (Fig. 2]) as peak 3, which corresponds to the decay 
of the level to the first excited state. It is predicted to de- 
cay by proton emission as well, but due to contamination 
from f3 particles and electronic noise it is not possible to 
identify it in Fig. [3] /3 particles mainly deposit a small 
amount of energy in the DSSSD and give a larger sig- 
nal in the pad detector. Since this source of background 
cannot be identified in the DSSSDs without backing, they 
are omitted in the following analysis. 

The large background below 5 MeV in Fig. [3] consists of 
multiplicity two events where one of the signals is caused 
by a /? or electronic noise typically in coincidence with a 
proton from a strong one- proton peak, e.g., a proton from 
31 C1 to the ground state in 30 S (see Fig. [J). The back- 
ground can thus be reduced by including only multiplicity 
two events with a proton that is known to feed the level 
considered. The protons feeding the 4689-keV level are 
found by gating on the 2478.9-keV 7. The protons are 
found to have energies in the intervals [1580, 1920]keV 
and [2200, 2400]keV. By gating on protons with these 
energies the background is reduced substantially. The 
4689-keV level, however, is still not positively identi- 
fied in the two-proton spectrum and thus only an upper 
limit on the r p /r 7 ratio can be found. The background 
can be estimated by gating on protons in the interval 
[2040, 2120]keV; they correspond to a known strong pro- 
ton group at 2084(2) keV from 31 C1 to the ground state 
of 30 S [ljl. The background spectrum then has to be 
scaled by a factor / = 0.311(4), namely to the ratio be- 
tween the number of protons gated on in the interesting 
intervals compared to the number of protons gated on 
for the background. The number of two-proton events, 
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TABLE III. Ratios between calculated proton and 7 partial widths from previous work compared with the 95 % confidence 
limits extracted from the present work. 



Setoodehnia et al. [10] 


A 1 

Almaraz- 


-Calderon et al. [14] 


Present work 


E x (keV) 






E x (keV) 


J 7 * 


r 7 /r p 


e x (keV) r 7 /r p 


4699(6) 


3+ 


210 








4689.2(24) > 3.8 


4814(3) 


2+ 


1.3 










5136(2) 


3+ 


480 


5130 


4+ 


340 




5217.8(14) 


0+ 


> 0.36 x 10" 3 


5217.8 


0+ 


0.7014 x I0~ 3 


5227(3) < 0.5 


5314(7) 


3- 


9.8 x 10~ 3 


5312.1 


3- 


4.6 




5391(3) 


2+ 


0.26 x I0~ 2 


5382 


2+ 


3.8954 x 10~ 2 










5835.5 


4+ 


15.7 


5847(4) < 9 



giving a 30 S energy in a 40-keV interval around 4689 keV, 
is found to be (26 — /33). Using the same proton gates, 
the number of 7 rays, corresponding to the decay of the 
4689- keV level to the first excited state in 30 S, is found 
to be (13-/11). 

Due to the limited number of counts in the spectra 
it is necessary to determine the upper limit through a 
simulation. The upper limit is found from an ensemble, 
which is made by drawing four numbers rij from a proba- 
bility function P(ni\noi), where noi ar e the four numbers 
found from the spectra (noi = 26, hq2 = 33, 7J03 = 13, 
rio4 = 11). This probability function is derived as fol- 
lows: n oi is a number from a Poisson distribution with 
mean value A^ and the n^'s needed to create the ensemble 
are random numbers drawn from a Poisson distribution 
with this mean value Ai. The problem arises because Ai 
is unknown. Therefore, one has to integrate over all the 
possible values of Ai weighted with the probability of this 
Aj given n 0i : 

P(ni\n 0l ) = oo 

J dX i P{X l \n Ql ) 

= / dAiPfalAiJPfnoilAO (2) 
Jo 

_ (nj+noi)! 1 

~ riilnoil 2 n *+ n °* +1 ' 

where the second equality comes from Bayes's theo- 
rem that employing a uniform prior distribution gives 
P(Xi\rii) = P(ni\Xi) and from the Poisson distribution 
being normalized in Ai. 

The r p /r 7 ratio can then be found from these four 
random numbers m drawn from the four distributions 
by correcting for efficiencies: 

-Tp = (ni ~ fn 2 )/e p ^. 
T 7 ("3 - fni)/e 1 ' 

where the 7 efficiency e 7 includes the 7 intensities mea- 
sured by Lotay et al. [13] . From the ensemble of 10 6 
r p /P 7 values found in this way, we find a 95 % confi- 
dence upper limit of 0.26 on the T p /T 1 ratio. 

This limit can be compared to calculations made by 



Setoodehnia et al. [10| for a 3 + resonance at 4699 keV. 
They find: r p /r 7 = 0.47 x 10~ 2 . Recent developments 
with nano structured CaO targets have increased the 
yield of 31 Ar up to an order of magnitude [25j]. Hence, 
using this type of target and a setup optimized for low- 
energy protons, it should be possible to identify the pro- 
ton decay of the 4689-keV level using the gating tech- 
nique presented here and thereby deduce an experimental 
value for the T p /T 1 ratio. 



The 5227-keV level 

A level at 5218.8(3) keV has been observed to 7 decay 
to the first excited state of 30 S by a 3008.5(2)-keV 7 and 
to the second excited state by a 1814.4(3)-keV 7 with 
relative branching ratios of 0.80(9) and 0.20(12), respec- 
tively [l3j. These lines are not seen in the 7 spectrum. 
Assuming this is the same level as we have identified in 
the two-proton spectrum at 5227(3) keV, an upper limit 
on the r 7 /r p ratio can be found by gating on the pro- 
tons feeding the level. These protons are found by gating 
on the 5227-keV peak in the two-proton spectrum and 
choosing the proton with the highest energy as this is 
most likely to be the first emitted proton. The two de- 
tectors without backing are again excluded. The protons 
feeding the level are found to have energies in the inter- 
vals [1000, 1350] and [1670, 1950]keV. The background is 
estimated by gating on another strong transition to the 
ground state of 30 S with proton energies in the interval 
[1380, 1460]keV (the proton group at 1416(2) keV in Ref. 
!16j). The number of two-proton events in the 30 S peak 
at 5227 keV is 226 ± 15. The are no 7 rays in a 50-keV 
interval around 3008 keV, corresponding to the 7 decay 
to the first excited state in 30 S. This gives a 95 % upper 
value for the number of 7 rays of 2.7 [26|. Including the 
7 intensities from Ref. [13[ , this gives a 95 % confidence 
upper limit of the r 7 /r p ratio of 0.5. 

As seen in Table UTT] Setoodehnia et al. [TtJ and 
Almaraz-Calderon et al. [l4[ estimated for a + level 
at 5217.8 keV that r 7 /P p > 0.36 x 10" 3 and r 7 /r p = 
0.7014 x I0~ 3 , respectively. Almaraz-Calderon et al. [l4| 
found a proton branching ratio of 1.00(2) for a level at 
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5200(44) keV, while Lotay et al. 13[ observed 7 decays 
of a level at 5218.8(3) keV. These results are hard to 
reconcile if the same state was populated. It should be 
noted that the calculations by Setoodehnia et al. and 
Almaraz-Calderon et al. are made for a + state, but it 
is not clear whether the level observed around 5220 keV 
is + or 3 + , or if there are in fact two levels in this energy 
region. 



The 5847-keV level 

The level at 5847(4) keV has been observed to 7 decay 
to the hrst excited state of 30 S by a 3637.7(4)-keV 7 ray 

.13j- There is no clear evidence of such a line in the 7 
spectrum and an upper limit on the r 7 /T p ratio can thus 
be found as described for the 5227-keV level. The two de- 
tectors without backing are again excluded. The protons 
feeding the level are found to have energies in the inter- 
vals [2220,3400] and [5540, 6160]keV. The background is 
chosen as for the 4689-keV level. The number of two- 
proton events in the 30 S peak at 5847 keV is 24 ± 5. The 
number of 7 rays in a 50-keV interval around 3638 keV 
is (2 - /12), where / = 0.463(4) defined as above. This 
gives a 95 % upper value for the number of 7 rays of 2.88 

26]. After correcting for efficiencies this gives a 95% 
confidence upper limit of the T 7 /r p ratio of 9. 

Almaraz-Calderon et al. [l4| do not observe any signifi- 
cant proton branch from this level in their 28 Si( 3 He, n) 30 S 
experiment. This does not agree with our result of a pro- 
ton branching ratio of at least 0.1 (95% C.L.). They as- 
sume it to be a 4 + state and estimate that: r 7 /r p = 15.7, 



in contradiction to our upper value of 9 (95% C.L.) 



IV. SUMMARY 

The levels below 6MeV in 30 S have been studied by 
the /32p- and /3p7-decay of 31 Ar. 

The 7 decay of the second excited state of 30 S to the 
ground state has been observed for the first time in the 
decay of 31 Ar. The relative intensities of the 7 lines corre- 
sponding to this decay and the decay to the first excited 
state have been found to be 0.63(21). 

The 7 decay of the astrophysically interesting level at 
4689.2(24) keV has also been observed for the first time 
in the decay of 31 Ar. 

We present a new analysis method that provides exper- 
imental limits on the ratio between the proton and 7 par- 
tial widths of resonances in 30 S. The upper limit of the 
T p /r 7 ratio has been found for the level at 4689.2(24) keV 
to be 0.26 (95 % C.L.), and upper limits of the r 7 /r p ra- 
tio of 0.5 and 9 (95% C.L.) have been found for levels 
at 5227(3) keV and 5847(4) keV, respectively. The latter 
conflicts with previous calculations. 
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